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ABSTRACT 

The axisymmetric linearized Euler equations are used to simulate noise amplification and radiation from 
a supersonic jet. The effect of heating on the noise field of the jet is studied and compared to experimental 
results. Special attention was given to boundary treatment, and the resulting solution is stable and nearly 
free from boundary reflections. 

NOMENCLATURE 


D 

jet exit diameter 

E 

mean total energy 

F, G 

inviscid flux vectors 

M 

Mach number 

P 

mean pressure 

Q 

vector of conserved properties 

T 

mean temperature 
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Subscripts 


U 

V 

Voo 

b 

e 

h 

P' 

u' 

v' 

t 

x, r 

a 

At 

Ax, Ar 

7 

co 

P 

P’ 


mean axial velocity 
mean radial velocity 

uniform mean radial velocity in the outer region 

half-width of the annular mixing layer 

perturbation total energy 

radius of the uniform core 

perturbation pressure 

perturbation axial velocity 

perturbation radial velocity 

time 

axial and radial coordinates 

complex wavenumber of the inflow disturbance 

time step 

grid spacing in the axial and radial directions 

ratio of specific heats 

frequency of the inflow disturbance. 

mean density 

perturbation density 


c 

j 

o 

t 

X, r 


oo 


centerline value 
jet exit centerline value 
stagnation value 
temporal operator 
axial and radial operators 
value at infinity 
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Superscripts 


spatial index 
time level index 


i 

n 


perturbation value 


INTRODUCTION 

The full, compressible Navier-Stokes equations govern the process of sound generation and propagation 
to the far field. To solve these equations without resorting to modeling the turbulent quantities. Direct 
Numerical Simulation (DNS) may be employed. However, current computer limitations make DNS 
calculations impractical, due to the large number of computational points required to resolve the high 
Reynolds number turbulent flows found in the shear layer of a supersonic jet. 

Therefore, Mankbadi et. al . 1 - 2 explored the use of the large-eddy simulation (LES) approach for the 
prediction of sound generation and propagation. In this approach, the full Navier-Stokes equations are 
divided into large-scale components, which are calculated directly, and small-scale components, which are 
unresolved, and thus modeled. The LES approach does not capture the sound radiation from the 
unresolved scales; however, it is believed that the large scales are much more efficient than smaller ones 
for radiating sound. Thus, LES is currently the most accurate approach to jet noise predictions that is 
computationally feasible. However, the LES approach is still CPU-intensive, especially for three- 
dimensional computations of both the near and far field. 

In previous work, the authors explored the use of the much less computationally-demanding linearized 
Euler equations (LEE) for jet noise predictions . 3 The LEE approach neglects both viscosity and nonlinear 
effects. The viscous effects can be safely neglected since the large-scale dynamics in free shear flows are 
essentially inviscid . 4 Nonlinearity, however, seems to be important . 5 
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Much of the physics of the problem are contained in the linearized equations. Several attempts have 
succeeded in studying the physics of jet noise based on a simplified form of the linearized Euler 
equations. 6 ' 9 Since the linearized Euler equations solve for the near and far field flow simultaneously, the 
problem of matching the sound generation in the near field to the propagation of sound to the far field does 
not arise. Also, the linearized Euler equations fully account for non-parallel flow effects and for the 
simultaneous presence of non-discrete frequencies. 

In the current work, the problem of sound generation and propagation from a heated jet is investigated. In 
cold supersonic jets, the structure of turbulence resembles a wave-like pattern (e.g., Ref. 1, 10). Thus, it 
is believed that the classical Kelvin-Helmholtz instability wave can model the large-scale structure in a 
supersonic jet, and the radiated sound is calculated accordingly (e.g.. Ref. 7). However, for heated jets at 
high Mach numbers the experimental results of Oertel 11 have shown that the jet can sustain additional 
instability modes in addition to the regular Kelvin-Helmholtz mode. 

The effect of these additional modes on the radiated sound has been addressed by Tam et. al. 12 (1992) and 
Seiner et. al. 13 (1993) using the Rayleigh equation. The Rayleigh equation can be obtained from the 
linearized Euler equation by assuming a locally-parallel flow and that the solution takes the form of normal 
decomposition. The radiated sound can then be determined via an asymptotic matching approximation of 
the instability wave solution to the acoustic field, as in Ref. 7. 

In the present work, the full linearized Euler equations are solved, which removes the need for the above 
assumptions. It is expected that the linearized Euler equations can capture the effects of all three families 
of instability waves if they are present. 

GOVERNING EQUATIONS 

Starting from the full Navier-Stokes equations in conservative form, neglecting viscosity, and linearizing 
about a mean flow (U,V), the axisymmetric linearized Euler equations may be written in cylindrical 
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coordinates as: 


where: 


dQ dF 1 d(rG) 1 a 
dt dx r dr r 


( 1 ) 


P 

6 = 1 -i, 

V 

e 

u 

p' +20U - pU 2 
uV + vU-pUV 
(p'+e)U + (u-pU)E 


G = 


v 

OV +vU -pUV 
p'+2vV-pV 2 
(p'+e)V + (v - pV)E 


( 2 ) 


( 3 ) 


( 4 ) 


and 


Here 



( 5 ) 


e - (OU + vV) 


p'=(y- 



p(u 2 + V 2 ) 


p = (/-1)p 



(6) 

( 7 ) 
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and 


(p,u,v,e) = [p' ,{pu)' ,{pv)' ,(pe)']. (8) 

Velocities are normalized by the jet exit velocity, time by D/Uj, density by the mean exit value at the jet 
centerline, and pressure by pjU_j 2 . 

MEAN FLOW 


To apply the linearized Euler equations, the mean flow field of the jet must be specified. In Seiner, et. 
al. 14 , the jet flow is expressed by means of a half-Gaussian profile described by the parameters b(x), h(x), 
and U c (x). 


There are three regions in the jet exhaust. The first two are the potential core and the transitional region, 
where the mean flow profile is described by: 


U = U c (x ) 


U = l/ c (x)exp 


-ln(2)( 


V 


r-h(x) 

b(x) 


>2 


J 


r < h 
r > h 


(9) 


In the potential core, U c (x) is equal to Uj, while in the transitional region, the centerline velocity is 
decreasing. 


At the end of the transitional regime, the parameter h(x) becomes zero, and the fully-developed regime is 
realized in which: 

U = U c {x)exp 


-ln(2) 


' r 


b{x) 


( 10 ) 


6 



The parameters h(x) and U c (x) are related to b(x) through the conservation of axial momentum: 

°° 1 

\pU 2 rdr = ^pjUjD 2 ( 11 ) 

0 


Using the boundary layer assumption, the pressure across the jet is taken as constant. The density in the 
jet is related to the mean velocity using Crocco's relation: 


EL = 
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( 12 ) 


In the potential core, U c = 1.0, and Eq. (1 1) is used to obtain h(x) in terms of a given b(x). For the fully- 
developed region, h(x) = 0.0 and Eq. (11) is used to obtain U c (x) in terms of a given b(x). In the 
transitional regime, a cubic spline fit is used to match the values of h(x) and its derivatives, then Eq. (11) is 
used to obtain U c (x) in terms of a given b(x). 

Once the axial velocity and density are known, the continuity equation for the mean flow can be integrated 
to obtain the mean radial velocity distribution. 

The above profiles are used to describe the mean flow up to a maximum radius r max = h + 3b. For r > 

^max- 

u = 0 

v < 13) 

v = — 

r 

The computational grid for this problem extends axially from x/D = 2.5 to x/D = 35, using 196 equally 
spaced points (approximately 25 points per wavelength). The computational grid starts 2.5 diameters 
downstream of the jet exit due to the numerical problems associated with the steep mean-flow gradients 
near the nozzle. 
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In the radial direction, the grid begins just above the centerline (r/D = 0.01) and extends to r/D = 16, with a 
total of 381 points. The grid is uniform from the centerline to r/D = 1, with a spacing of Ar/D = 0.01. At 
this point, the grid is stretched geometrically by a factor of 1.01 until the radial spacing is equal to the axial 
spacing. After this point, the grid is uniform again to the outer radial boundary. 

NUMERICAL FORMULATION 


The code is a modified split MacCormack-type solver, which is second order accurate in time and fourth 
order accurate in space. This extension of the MacCormack scheme is known as the 2-4 scheme, and was 
developed by Gottlieb and Turkel. 15 This scheme has been used successfully on a wide range of unsteady 
fluid and aeroacoustics problems (e.g. Refs. 2, 3, 8, 16-22). Sankar et. al 23 have evaluated this scheme 
for aeroacoustics applications. This code has previously been validated on an unheated jet problem by 
Mankbadi et. al. 3 . The solution procedure is as follows: 

In the present code, the operator is split into separate radial and axial contributions: 

Q n+2 = L x L r L r L x Q n (14) 


Each operator consists of a predictor and a corrector step. Each step uses one-sided spatial differencing: 


Predictor: 


1 

- n+ ~ 

Q 2 =Q n 



Corrector: 


6 n+1 


2 

2 


1 

- n+ “ 
Q n + Q 2 


l + 


At 
6 Ax 


(7F,-8F; + i 


A 



(15) 


(16) 


8 



The radial contribution is computed in a similar way. The order of the sweep directions are reversed 
between operators to avoid biasing. At the computational boundaries, flux quantities from outside the 
boundaries are needed to compute the spatial derivatives, and these can be obtained by using third-order 
accurate extrapolation based on information from the interior of the domain. 

The boundary treatment used is similar to previous work. 3 On the inflow boundary in the hydrodynamic 
disturbance regime (x/D < 2), a characteristic-based boundary condition is specified in the jet flow. 24 - 25 
In the still air of the outer radial boundaries, an acoustic radiation treatment is specified. On the jet outflow 
boundary, the asymptotic outflow condition of Tam and Webb 26 is used. At the centerline, an averaging 
procedure is used. 27 A more complete description of the implementation and the boundary conditions 
used may be found in Refs. 2 and 28. With these boundary conditions, the solution is stable and free of 
reflection from the boundaries. 


At the inflow boundary of the computational domain, a small disturbance is introduced. This disturbance 
is assumed to be mainly hydrodynamic in nature, and is specified from the centerline to r/D = 2.5. 


To a first approximation, the inflow disturbance is assumed to be small such that the linear stability theory 
applies. A normal mode decomposition for the disturbance is assumed to be of the form: 


lu'y,p\p'\= 

Re[t/(r ),v(r ),p(r ),p(r )]exp{/(ax - art)} 


(17) 


The governing equations reduce to the Orr-Sommerfeld equation, which is solved to obtain the complex 
wave number oc as the eigenvalue corresponding to the frequency (0 and the radial functions (~) as the 
corresponding eigenfunctions. The mean flow discussed above is used in solving the Orr-Sommerfeld 
equation. This solution extends to r/D = 2.5. A curve is fitted to smoothly reduce the disturbance to zero 
by r/D = 3. 
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The effect of the inflow disturbance is added to the computed flow variables at the inflow boundary at each 
time step: 



boundary 



computed 


+ 



disturbance 


RESULTS 


08) 


Results are presented for the axisymmetric flow and acoustic field of unheated and heated supersonic (Mj 
= 2.0) jets, with a farfield temperature of 279.5° K. Three cases were computed: (1) a T„ = 313° K jet 
excited at a Strouhal number of 0. 19, (2) a T 0 = 755° K jet excited at a Strouhal number of 0.22, and (3) a 
T 0 = 1 1 14° K jet excited at a Strouhal number of 0.2. These cases were tested experimentally by Seiner, 
et. al. 14 

The axisymmetric code was run on the NAS Cray C90, and required 4.5 mW of memory and 30 minutes 
of CPU time for each run. 

Figures 1, 2, and 3 show the instantaneous distribution of pressure at a time level during the computation. 
In all three figures, the dominant angle of radiation is clearly illustrated. As the core temperature of the jet 
is increased, the angle of emission becomes steeper. Also, as the temperature rises, the supersonic family 
of waves appear, and sound is radiated at two dominant angles. Figures 2 and 3 show an interference 
pattern due to the appearance of this second family of instability. 

Figures 4, 5, and 6 show the root-mean-square pressure levels. Again, it is seen that the noise emission 
angle steepens with increasing core temperature. Also, the interference pattern due to the appearance of the 
supersonic instability is seen in Figs. 5 and 6. 

Figures 7, 8, and 9 compare the directivity of the sound field to the experimental results of Seiner, et. al. 14 
It should be noted that the distance of the computed results from the nozzle is much smaller than that of the 
experiment. The experimental data was gathered between 40 < R/D <72.3, while the computed data was 
gathered between 16 < R/D < 36.25 (due to the limits of the computational domain). This will cause the 
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computed radiation angles to be approximately 3° lower than the experimental values. Since the amplitude 
of the disturbance at the jet exit plane is not known, the peak of the sound pressure levels are scaled to 
match that of the experimental data in each case. In these figures, a peak is shown at the 100 degree angle, 
the first three figures show that this peak is an artifact of the input disturbance, which travels up the inflow 
boundary. 

In these figures, it is shown that the computation captures the angle of maximum noise radiation well, even 
though only the axisymmetric (n=0) mode is being computed. The steepening of the angle of maximum 
noise radiation with increasing temperature is also correctly predicted. 

CONCLUSIONS 

A linearized Euler equation approach is used to study the effect of increasing core temperature on the 
sound propagation in a supersonic jet. The axisymmetric (n=0) mode is computed, and the results are 
compared with experiment. Good qualitative results are obtained for the angle of maximum noise radiation 
and for the shift in the angle of maximum noise radiation as the jet temperature is increased. 
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Figure 2. — Instantaneous pressure distribution (T 0 = 755 °K; St = 0.22). 
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Figure 3. — Instantaneous pressure distribution (T 0 = 1114 °K; St = 0.20). 



Figure 4.— Sound pressure level distribution (T 0 = 313 °K; St = 0.19). 
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Figure 5. — Sound pressure level distribution (T 0 = 755 °K; St = 0.22). 



Figure 6. — Sound pressure level distribution (T 0 = 1114 °K; St = 0.20). 
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Figure 9. — Sound pressure level comparison (T 0 = 1114 °K; St = 0.20). 
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